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Resistive leakage current based condition assessment of metal oxide surge 
arrester (MOSA) is one of the most extensively employed technique to 
monitor its degradation. An extraction method is customarily required to 
extract the resistive component from the total leakage current. The existing 
methods to extract the resistive current are complex and less accurate. 
Therefore, this paper describes a simple and accurate circuit-based method to 
extract the resistive current using equivalent model and measured leakage 
current of the arrester. The accuracy of the proposed method is validated 
through experimental results on ABB’s 120 kV surge arrester, EMTP 
and QuickField software simulations. The performance of the method is also 
analyzed and verified experimentally on 72, 180 and 240 kV rated ABB’s 
surge arresters. The obtained results of resistive leakage current have shown 
the maximum error of 0.001%. Simple and easier computational steps with 
higher accuracy are the key benefits of the proposed technique. 
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1. INTRODUCTION 

Surge arresters are employed to protect the transmission system from the over voltage caused 
by lightning and switching operations in the power system [1]. The general types of arresters are silicon 
carbide-based arresters with spark gaps and gapless metal oxide arresters. Gapped arresters are not able to 
limit the switching overvoltage effectively [2]. Therefore, the gapless arresters are developed; the most 
popular MOSA is the zinc oxide arrester [3]. Zinc oxide element offers high nonlinear voltage current 
characteristics and faster conduction response for high voltage surges [4, 5]. As the resistance of zinc oxide 
element cannot be made infinity, therefore it draws a continuous current at the normal operating voltage, 
known as leakage current [6, 7]. The arrester’s total leakage current consists of the resistive and capacitive 
components [8]. The resistive current varies with the changes in zinc oxide’s characteristics, operating 
voltage, ambient temperature, and environmental factors [9-12]. Therefore, it can be used as the most reliable 
indicator of ageing and degradation of MOSA [13]. 

Many techniques have been proposed for the extraction of resistive leakage current which include; 
capacitive current compensation method [14, 15], field probe method [16, 17] and current orthogonality 
method [4, 18-21]. In the current compensation method, the capacitive current is compensated by adjusting 
the amplitude and phase angle of the reference signal but the presence of harmonics in the capacitive current 
may interfere with the measurement of the nonlinear resistive current. In the field probe method, the third 
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harmonic capacitive current by voltage harmonics is compensated with the help of a compensating 
current signal but the effect of voltage distortion on third harmonic resistive current is not considered. 
Similarly, the modified shifted current method [18-20], time delay addition method [21], and current 
orthogonality method [22] are proposed, which are based on the displacement of the resistive and capacitive 
components of the leakage current. Modified shifted current and time delay addition methods are based on 
waveform manipulation of the obtained leakage current. A capacitive current waveform with a shifted phase 
of a radians is injected to the total leakage current to compensate the capacitive current component. 
The magnitude of the capacitive current is determined at the time instant, which has a time delay of 1/2 
radians with respect to the peak time of resistive current waveform [23]. However, the results of acquired 
resistive current are more accurate as compared to the conventional compensation technique. 
But, the accuracy of results is affected in the case of non-sinusoidal applied voltage [24]. Lot of complex 
and complicated steps are required to extract the resistive current in these methods and the results are not 
enough accurate [25]. Therefore, the implementation of an extraction method with simpler computational 
steps and better accuracy provides a potential research area. The aim of this paper is to propose a simpler 
and accurate computational circuit-based method for the extraction of the resistive leakage current of MOSA. 
The proposed method is based on the measurement of total leakage current and mathematical equations 
derived from the equivalent model of MOSA. The accuracy of method is validated experimentally on ABB’s 
120kV surge arrester and through particulars software such as EMTP and QuickField. Its performance is also 
tested experimentally on 72, 180 and 240kV surge arresters. The obtained results have revealed that the 
proposed method has higher accuracy with minimum error of 0.001%. Furthermore, the implementation 
procedure is quite easy and simpler as compared to other existing methods. 


2. RESEARCH METHOD 
2.1. Equivalent circuit of MOSA 

An equivalent circuit model of MOSA as shown in Figure 1, is required to determine the circuit 
element values and resistive leakage current [26]. In Figure | ‘V,’, ‘I? and ‘I, represent the supply voltage, 
total and resistive leakage current respectively. Equestion (1) presents the orthogonal relation between 
the resistive and capacitive components of the total leakage current. 





Figure 1. An Equivalent circuit model of MOSA 


LELEI (1) 


The capacitance ‘C? and resistance ‘R’ in the circuit model depend on the electrical properties and 
dimensions of the MOSA column. ‘C’ can be determined using (2). 


E9ErA 
d 


C= (2) 


Where; 
£» = Vacuum permittivity 
& = Relative permittivity of the arrester disc 
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A = Area of the MOSA disc 

d = Height of the column of MOSA disc 

Furthermore, the resistivity ‘p’ and resistance ‘R’ of the MOSA disc can be represented as a function of 
the applied voltage ‘V,’ in (3) and (4). 


p=f(%) (3) 
RWVs) = 22" (4) 


2.2. V,- I characteristic curve of the surge arrester 

The voltage-current characteristic curve, as shown in Figure 2, shows the non-linear variation of 
applied voltage and the total leakage current. Three regions are defined in VY, — J, curve based on applied 
voltage. Four points are destined in the region I, known as the steady state region. At the first point (Vc), 
the total leakage current of the MOSA is mainly capacitive. While the other points are the phase to earth 
voltage (V,_y), continuous operating voltage (Vcoy) and rated voltage (Vrateq). A small increase in 
the voltage increases the leakage current significantly in the region II, called the flat region. The temporary 
over voltage (Vrovy) also occurs in region IT. Whereas, in the region III known as an arrester operation region, 
MOSA operates to discharge the lightning current. The discharge voltage, which appears in this region 
is called ‘Vj’. The considered values of the voltages such as V., Vp.n, Vcov and Vyateq are presented in 
the Table 1. 


Table 1. Considered value of V, ,Vp_y, Veoy and V;ateq for 120kV surge arrester 








Parameter Value 
V; TOK Vrms 
ae 93k Vis 
Vtg 98k Vis 
V,-atea 120 kVims 





Voltage (V,) 





Leakage current (I) 


Figure 2. V, — I, characteristic trend of a typical MOSA 


2.3. Extraction of the resistive leakage current using mathematical procedure 

The relationship between the applied voltage and total leakage current based on the equivalent 
circuit of Figure 1 is given in (5), where ‘Z’ is the impedance of the equivalent circuit of MOSA. 
At the voltage levels, V, < V,, the contribution of the resistive component in the total leakage current 
is negligible. Therefore, the impedance of the circuit model is given in (6). 


V =Z (5) 


1 


a ae (6) 


Where; 
X= Capacitive reactance of the circuit 
œ= Angular Frequency 
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The capacitance of the MOSA disc (shown in 7) can be determined by substituting (5) and (6) 
into (1) and by considering ‘I,’ as ‘I,’ and ‘V,’ as ‘V,’. Similarly, the capacitive leakage current in the steady 
state condition (as shown in (8)) is obtained by equating (2) and (7). The resistive leakage current can easily 
be calculated after determining the capacitive current. 


Ic 





C= (7) 
I, = = (V, x 2r x 50) (8) 


2.4. Experimental procedure 

The experimental setup to measure the total leakage current of ABB’s 120kV surge arrester 
is shown in Figure 3. The thickness of each disc of the arrester is 2.5 cm while the height of the column of ‘d’ 
and radius 'r’ of the disc are 100 cm and 3.5 cm respectively. A high voltage transformer (HVT) is connected 
in series with a pure resistor of 970 Q to apply voltage to the tested arrester. The voltage across the shunt 
resistor is used to determine the total leakage current l}. The magnitude and waveform of the peak voltage 
across the resistor is measured and recorded by digital oscilloscope (Osc). 









HV Transformer 
Resistive divider 


Resistor 


Figure 3. Schematic diagram of experimental setup 


2.5. Validation of proposed method by EMTP software 

In order to validate the results of the proposed technique, initially the circuit model of MOSA 
as shown in Figure | is designed in EMTP software. The capacitance of the equivalent model is calculated 
using (2). The non-linear resistance of the arrester is determined using the curve fitting technique based on 
the plotted V,-J, curve of Figure 4. The resistive leakage current as a function of applied voltage is then 
computed using (9). The obtained resistive leakage current using EMTP software is compared with 
the results of proposed method and is presented in the result’s section. 


_ Vs 

p=% (9) 
2.6. Validation of proposed method using quickfield software 

The electrical properties such as relative permittivity (€,), and resistivity (p) are required for finite 
element modelling in the software. The obtained resistive leakage current using QuickField software 
is compared with the results of proposed method and is presented in the result’s section. 


3. RESULTS AND DISCUSSION 

The proposed method is applied to compute the capacitive and resistive components based on 
the measurement of total leakage current in the experimental test described in section 2.3. The obtained 
results for a 120 kV MOSA are presented in Table 2. The data in the first three columns of Ving Vpeay and I, 
is obtained from the experimental test. While, the last two columns of J, and J, are computed using 
the proposed technique. By analysing the values of J, 7, and J,, it is found that the leakage current 
is dominantly capacitive for the voltage level up to the continuous operating voltage of 98 kVims. 
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However, the resistive current significantly increases for voltages above 120 kVins. The variation of J,, I, 
and J with the applied voltage can be seen in Figure 4 based on the values from Table 2. The Y%-I, 
characteristic is linear for the applied voltages in the range of 14.14 to 131.52 kV peak as shown in Figure 4. 


It can be said that Vp_y is the knee point of the Y.-J, curve of the arrester. After Vp_y, the resistive current 
shows a non-linear characteristic. Whereas, the variation of the capacitive current is linear for all ranges of 
the applied voltage. 
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Figure 4. The variation of resistive current (I,), capacitive current (I,) and total leakage current (I, ) 
of 120kV MOSA with the applied voltage using the proposed computational method 


Table 2. Computed capacitive and resistive leakage current of 120 kV surge arrester using the proposed 
mathematical method 








Vas Vpeak l i i 
(kV) (kV) (uA) (WA) (WA) 
10 14.14 149.48 149.48 0.80 
20 28.28 298.98 298.96 1.61 
30 42.42 448.56 448.45 2.41 
40 56.56 598.14 597.93 3.22 
50 70.71 747.94 747.41 4.02 
60 84.85 897.94 896.89 4.82 
70 98.99 1047.94 1046.38 5.63 
80 113.13 1197.94 1195.86 70.55 
93 131.52 1393.81 1390.19 100.51 
98 138.59 1474.23 1464.93 165.33 
110 155.56 1731.96 1644.31 544.00 
120 169.70 2680.41 1793.79 1991.72 
130 183.84 10124.74 1943.27 9936.50 





3.1. Validation of proposed method by EMTP software 

Based on the data presented in Table 3 to design the equivalent circuit model of MOSA, 
the capacitance is found to be 33 pF. The resistive leakage current as a function of the applied voltage at 
the steady state condition by using the curve fitting technique is given in (10), which can be substituted 
in (9) to determine the resistance. The voltage dependent function of resistance R(V,) is shown in (11). 
The obtained values of resistance at V, ,Vp_y, Vcoy and V-ateq are presented in Table 4. The peak values of 
the total leakage current and its resistive and capacitive components are presented in Table 5. 


1, (Vz) = 0.0038 x e29774%s (10) 
_ Vs 
RCV) ~~ 9.0038xe0:077405 (1) 


Table 3. Computed data for 120 kV surge arrester 











Parameter Formula Value 
C ExA 33pF 
d 
E (A 8.68 (nF/m) 
A nr? 0.003848 (mô) 
d - 1 (m) 
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Table 4. Resistance of surge arrester at different voltages 








Applied voltage (V, ) Applied voltage (Y, ) Resistance of the 
(rms value) Peak value arrester 
RV) 
10kV;ms 14.14 kV 1.24T 
93KV;ms 131.52 kV 1.31 GN 
98kVims 138.59 kV 800 MQ 
120 KVims 169.7 kV 88 MQ 





Table 5. Simulation results of 120 kV surge arrester modeled in EMTP 








Vems Vpeak I, I, L 
(kV) (kV) (WA) (WA) (WA) 
10 14.14 147.95 146.6 0.80 
93 131.52 1379.53 1363.49 100.20 
98 138.59 1460.03 1436.83 173.18 
120 169.70 2617.86 1759.36 1924.55 





3.2. Validation of proposed method by quickfield software 

The relative permittivity can be derived by substituting the dimensions of the arrester in (2) and 
non-linear resistivity formula (derived from (4)) is given in (12). The electrical conductivity (o) can be 
determined by taking the reciprocal of resistivity as shown in (13). The obtained values of resistivity 


and conductivity at different values of applied voltages is presented in Table 6. The simulation results of 
the leakage current values are shown in Table 7. 


mr? x Vs “J 
0.0038dxe”?7 Vs 


(p(V)) = 


d 


1 
P 


nr? xR(Vs) 


(12) 


(13) 


Table 6. Electrical resistivity values of 120 kV surge arrester at different values of applied voltage 


Applied peak voltage (V; ) 


Electrical resistivity of the 


arrester p (GN. m) 


Electrical conductivity of 


the arrester o (s/m) 





14.14 kV 4.79 2.086x10"” 
131.52kV 0.005 1.9793x107 
138.59 kV 0.003 3.2469x107 
169.7kV 0.00033 2.9468x10° 





Table 7. Simulation results of 120 kV ZnO arrester modeled in quickfield 








Vis Vpeak I, I, L 

(kV) (kV) (WA) (WA) (WA) 
10 14.14 147.84 147.84 0.1 
93 131.52 1378.42 1374.72 100.99 
98 138.59 1458.23 1447.68 175.10 
120 169.70 2628.28 1774.08 1939.20 





3.3. Comparison between the computed and simulated results 
The implementation of the proposed technique yields computed results whereas the simulated 


results can be obtained by using software. Table 8 shows less than 1% difference between the obtained 
leakage current by computation and simulation using EMTP software. Moreover, the differences between 
the simulated and computed resistive and capacitive components of the leakage current are less than 5% 
and 2% respectively. This difference arises due to curve fitting technique to evaluate the resistive leakage 
current at continuous operating voltage. Similarly, Table 9 shows less than 1% difference between 
the computed and simulated resistive and capacitive leakage currents modelled in QuickField software. 
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4. 


Table 8. Comparison between simulated and calculated results of 120 kV surge arrester modeled 
in EMTP software 
L I 





I 








Vims (HA) (HA) (HA) 
Vv) Simulation Computation Simulation Computation Simulation Computation 
10 147.95 149.48 146.6 149.48 0.80 0.80 

93 1379.53 1393.81 1363.49 1390.19 100.20 100.51 
98 1460.03 1474.23 1436.83 1464.93 173.18 165.33 
120 2617.86 2680.41 1759.36 1793.79 1924.55 1991.72 





Table 9. Comparison between simulated and calculated results of 120 kV surge arrester modeled 
in QuickField software 











I, I, I, 

Vims (uA) (uA) (uA) 

(kV) Simulation Computation Simulation Computation Simulation Computation 
10 147.84 149.48 147.84 149.48 0.1 0.80 
93 1378.42 1393.81 1374.72 1390.19 100.99 100.51 
98 1458.23 1474.23 1447.68 1464.93 175.10 165.33 
120 2628.28 2680.41 1774.08 1793.79 1939.20 1991.72 





CONCLUSION 
In this paper, a circuit-based method for the extraction of the resistive leakage current has been 


proposed. The accuracy of the proposed computational method has been validated using 120 kV surge 
arrester’s experimental results and software such as EMTP and QuickField. Also, the performance of the 
proposed method has been compared and tested for ABB’s 72, 180 and 240 kV rated surge arresters. The 
obtained results have shown that the implementation of this proposed technique requires less computational 
steps with higher accuracy as compared to existing methods. The error between the simulated and calculated 
resistive current values has been found to be 0.001%, which is quite less than the errors obtained from the 
available methods. In the future, it is aimed to modify the proposed technique to extract the online resistive 
leakage current of MOSA. 
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